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SECTION  I 

INTRODUCTION  AND  SUMMARY 


Much  of  the  analysis  of  the  mantle  P-wave  noise  proposed 
for  this  contract  depended  on  the  ability  to  treat  the  LASA  as  a  large  array, 
i.  e.  ,  to  perform  coherent  processing  by  treating  seismic  noise  as  a  com¬ 
pilation  of  plane  waves  using  the  large  aperture  available  at  LASA.  The  sub¬ 
array  outputs  should  be  P-wave-limited,  and  the  large  array  could  be  used 
to  dissect  the  P-wave  noise;  therefore,  it  is  extremely  important  to  under¬ 
stand  the  coherence  of  the  subarray  outputs. 

Considerable  effort  has  been  expended  to  determine  coherence 
among  subarray  outputs.  This  report  presents  the  results  of  coherence 
studies  among  subarrays. 

The  principal  results  of  the  study  are 

•  Development  of  a  program  to  compute 
multiple  coherences  and  thus  estimate 
multichannel  coherence 

•  Development  and  programing  of  a  new 
processing  technique  to  design  filters 
that  maximize  the  coherence  between 
the  output  of  two  groups  of  sensors 
(group  coherence) 

•  Presentation  of  results  of  measured 
2-channel  and  multiple  coherence, 
indicating  that  the  subarray  outputs 
generally  are  moderately  coherent 
below  about  1  cps  only  within  the  B 
ring  of  the  subarrays 


/ 

/ 
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•  Presentation  of  2-channel  and  multiple 
coherences  for  several  possible  noise 
models.  A  study  of  these  coherences 
suggests  that  the  most  reasonable  models 
would  give  little  coherence  between  AO 
and  the  E  and  F  rings.  Comparison  of 
model-study  results  with  measured  co¬ 
herences  suggests  that  a  considerable 
portion  (25  percent  or  more)  of  the  sub¬ 
array  output  power  is  not  interpretable 
as  seismic  energy 
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SECTION  II 

DISCUSSION  OF  COHERENCE 


Coherence  is  a  measure  of  the  predictability  of  the  primary 
channel  from  the  other  channels  as  a  function  of  frequency.  Consider  a 
least-mean-square  prediction  of  one  channel  fro.  i  n- 1  others  in  the  frequency 
domain.  The  correlation  matrix  formed  from  the  variables  is  a  crosspower 
matrix  which  can  be  partitioned  as  follows: 


§  n,, 

11  12 


^21  ^ZZ 


where 


$  is  the  autopower  of  the  trace  being  predicted 


is  the  row  vector  of  crosspowers  $  through  $ 

12  12  In 


n21  18  ni2 


Q  is  the  crosspower  matrix  formed  by  channels  2  through  n 

w  w 


-1 


The  solution  is  filter  weights  of  Q  ^  ^2  anc*  Pre<liction-error 
variance  (autopower)  of  $  ^  -  0^  ^22  ^  ^21  or 


»  11- 
911  1 


ni2°22  n21 


$ 
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Coherence,  as  used  in  this  report,  is  defined  as 


§ 

11 


n 


21 


For  the  2-channel  case,  this  reduced  to 


$ 

? 


12 

11 


$ 

§ 


21 

22 


Thus,  (1  -  coherence),  whether  2-channel  or  multiple  coherence,  can  be 
interpreted  as  the  fraction  of  the  power  in  the  primary  trace  which  cannot 
be  predicted  from  the  other  traces.  True  coherence  is  unchanged  by  the 
application  of  linear  filtering  to  any  or  all  traces. 

Coherence  is  a  random  variable  which  must  be  estimated. 
Estimation  of  coherence  is  determined  from  the  estimation  of  the  cross¬ 
power  matrix.  The  crosspower  matrix  is  estimated  by  some  scheme  from 
a  finite  data  sample  and  may  not  accurately  reflect  the  ensemble  statistics. 
The  result  may  be  an  overestimate  or  an  underestimate  of  coherence. 


Schemes  to  estimate  power  spectra  involve  smoothing  (in  the 
frequency  domain)  to  give  the  estimate  stability.  This  smoothing  introduces 
a  bias  into  the  crosspower  spectral  estimates,  making  them  systematically 
low.  Smoothing  does  not  bias  the  autopower  estimates  and  thus  makes  the 
expected  values  of  the  coherence  estimates  too  low.  This  can  be  a  severe 
problem  when  there  is  considerable  time  delay  among  channels. 
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For  example,  consider  two  channels  which  are  exactly  the 
same  except  that  one  is  delayed  by  t: 

fj(t)  =  f(t) 

f2(t)  =  f  (t  +  T) 

Assume  that  the  traces  are  sampled  at  an  interval  of  AT  and 
that  the  samples  have  a  total  length  NAT.  Let  t  =  KAT.  If  the  crosspower 
matrix  is  to  be  estimated  by  direct  transform  methods,  the  data  are  trans¬ 
formed  according  to  the  dictates  of  the  sampling  theorem  (Af  =  l/NAT)  and 
the  crosspower  estimate  is  of  the  form 


a+C-l 

■—  ^  Fl  *  (MAf)  F2(MAf) 
M=a 


where 


Af  =  1/NAT 


C 

FL(f) 

F2(f> 


number  of  frequency  points  in 
smoothing  interval 

discrete  Fourier  transform  of 
fL(t) 

discrete  Fourier  transform  of 

yt) 


The  crosspower  at  each  frequency  (f  =  MAf)  will  be 


exp 


M 

NAT 


KAT 


exp 


MK 
N  - 


II- 3 


science  services  division 


For  notational  compactness,  let  K/N  =  Q.  Therefore,  the 
estimate  of  the  crosspower  becomes 


a+C-1 

~~  exp  (2n  i  MQ) 

M=a 


exp(2rTiaQ)  [exp(2TTiCQ)  -  l] 
C  [exp  (2tt  iQ)  -  l] 


sin  tt  CQ 
C  sinnQ 


exp 


2 


2tt  i  Q  (a  +  -y 


Thus,  the  crosspower  estimate  is  reduced  by  a  factor 
(sin  tt  CQ/C  sin  tt  Q). 

Consider  an  example  which  would  be  reasonable  for  two  LASA 
subarray3  (AO  and  C  ring)  with  100-point  smoothing.  For  a  noise  sample 
dominated  by  a  10-km/sec  plane  wave,  Kw  20,  N«s4000,  1/200.  and 

C=  100;  then,  the  coherence  would  be  reduced  by  about  1/1.  57. 

This  biasing  effect  makes  the  estimation  of  crosspowers 
between  the  widely  spaced  LASA  subarrays  very  difficult.  Extremely  long 
noise  samples  would  be  necessary  to  estimate  adequately  a  crosspower 
spectrum  from  two  positions  approximately  100-km  apart. 
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SECTION  III 

MULTIPLE-COHERENCE  PROGRAM 


A  rigorous  mathematical  derivation  of  the  multiple  coherence 
(R^)  for  weakly  stationary  Gaussian  processes  has  been  given  by  L.  H.  Koopmans 
The  calculation  is  based  on  the  fact  that  a  crosspower  matrix  is  a  covariance 
matrix  of  random  variables,  which  is  the  Fourier  transform  of  a  sample  of 
time-series  data.  Since  Fourier  transforms  are  linear  combinations  of 
multivariate  normal  data,  they  also  are  distributed  multivariate  normal.  One 
then  uses  the  body  of  theorems  about  conditional  distributions  for  multivariate 
normal  data.  The  actual  calculations  are  done  iteratively  (over  an  increasing 
number  of  sensors),  using  the  exact  equations  given  by  Koopmans. 

An  entirely  analogous  development  of  the  equations  for  ral- 

culating  multiple  correlations,  which  is  mathematically  much  easier  to 

2 

handle  rigorously,  is  given  by  T.  W.  Anderson.  This  method  of  estimating 
multichannel  coherence  or  predictability  is  computationally  very  rapid  com¬ 
pared  to  the  design  of  prediction  filters  and  the  evaluation  of  their  performance. 
The  method  is  especially  useful  when  sevexa.  ent  combinations  of  sensors 

are  desired. 

The  program  estimates  the  crosspower  spectra  by  taking  the 

direct,  transform  of  time-series  data,  using  the  Cooley-Tukey  algorithm, 

and  setting  «  (f)  =  F.  *  (f)  F.  (f).  Stability  is  obtained  by  averaging  adjacent 

ij  1  J  -  2 

frequencies.  The  distribution  of  the  estimate  of  the  multiple  coherence  (R  ) 

has  been  derived  by  N.  R.  Goodman3  under  the  assumption  that  the  frequency 

vectors  (direct  transforms  of  time-series  data)  used  in  the  averaging  are 

independent.  This  leads  to  density  function 


P(RZ) 


r(n) 

~  r (o-i) r(n-P+i) 


(1  -  R2)  (R2)P‘2  (1  -  R2)n'P  F  (n,  n;p-  1;R2  R  2) 
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where 


r 

F 

P 

n 


is  the  multiple  coherence  estimate 

is  the  true  population  multiple  coherence 

is  the  gamma  function 

is  the  hypergeometric  function 

is  the  number  of  channels 

is  the  number  of  frequencies  used  in  the  formation 
of  the  crosspower  matrix  estimate 


2 

This  reduces  to  a  Beta  distribution  when  the  data  are  uncorrelated  (R  =  0). 

In  estimating  the  multiple  coherences,  80-point  (0.  2-cps) 

smoothing  with  4096  (409.6-sec)  noise  samples  generally  were  used.  Table 

2 

HI- 1  gives  the  mean  for  R  =  0,  n  =  80,  and  p  =  1  through  15. 


Table  III  -  1 

MEAN  VALUE  OF  MULTIPLE- COHERENCE  ESTIMATE 


p 

Mean 

1 

1/82  =  0.  012195 

2 

2/82  =  0. 024390 

3 

3/82  =  0. 036585 

4 

4/82  =  0. 048780 

5 

5/82  =  0. 060976 

6 

6/82  =  0. 073171 

7 

7/82  =  0. 085366 

8 

8/82  =  0.  097561 

9 

9/82  =  0. 109756 

10 

10/82  =  0. 121951 

11 

11/82  =  0. 134146 

12 

12/82  =  0. 146341 

13 

13/82  =  0.  158537 

14 

14/82  =  0.  170732 

15 

15/82  =  0.  182927 
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0 

0 

D 

fj 

0 
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Burg  has  shown  that  the  ratio  of  the  mea.n-square-error  for 
an  estimated  filter  to  the  mean-square-error  for  the  true  ensemble  filter  in 
a  prediction  problem  is  independent  of  the  data's  covariance  matrix. 

—  2 

(1  -  R  )  can  be  interpreted  as  the  fraction  of  prediction- error 
power.  This  is  not  exactly  the  Urm  required  by  Burg's  theorem  (fractional 
as  opposed  to  true  prediction- error  power)  but,  if  the  data  sample  is  assumed 
to  give  an  accurate  spectral  estimate  of  the  primary  trace  (the  one  to  be  pre¬ 
dicted),  the  results  can  be  extrapolated  to  correlated  noise.  If  the  mean  value 
of  multiple  coherence  for  the  uncorrelated  data  were  0.  15,  then,  from 
[0.85  (MSE  in  prediction  estimated)/ 1.  0(MSE  for  ensemble  average)],  one 

would  expect  a  noise  sample  with  0.  98  true  multiple  coherence  to  give  an 
—  2  —  2 

estimate  (R  /0.  02)  «  0.  85.  Thus  1  -  R  s  (0.  02)  (0.  85)  a  0.  017. 

This  analysis  is  not  rigorous  but  does  indicate  that  the  esti¬ 
mated  multiple  coherence  becomes  much  more  reliable  (in  terms  of  absolute 
error)  as  the  true  coherence  increases. 
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SECTION  IV 


2-CHANNEL  COHERENCES 

To  get  an  initial  idea  of  the  coherence  among  subarrays, 
outputs  of  the  MCF  system  were  used  to  calculate  2-channel  coherences 
between  the  AO  and  the  Bl,  B2,  B3,  Cl,  C2,  C3,  C4,  Dl,  D2,  D3,  and 
D4  subarrays.  (Figure  IV - 1  shows  the  large-array  configuration.)  This 
was  done  for  an  8-min  noise  sample  recorded  25  March  1966  from  04:18  to 
04:26.  Data  were  prewhitened,  using  a  single  deconvolution  filter,  before 
the  coherences  shown  in  Figure  IV- 2  were  computed. 

These  coherences  are  probably  not  meaningful  at  very  low 
frequencies  (f  <  0.  1  cps),  because  nonseismic  noise  likely  makes  up  a 
considerable  portion  of  the  power  at  these  low  frequencies.  The  short- 
period  seismometer  response  is  such  that  low-frequency  ground  motion  is 
severely  attenuated. 

The  2-channel  coherences  have  the  following  salient  features: 

•  They  fall  off  sharply  with  increasing  distance; 
there  is  almost  no  coherence  between  AO  and 
subarrays  of  the  D  ring 

•  All  closer  subarrays  (within  13  km)  except  C4 
show  a  peak  coherence  near  1/3  cps,  which  is 
near  the  microseism  peak;  Figure  IV- 3  shows 
the  power  spectrum  of  the  MCF  output  for  sub¬ 
array  AO,  and  this  low-frequency  coherence  is 
very  likely  caused  by  the  directional  surface¬ 
mode  energy 

Theoretical  coherences  generated  by  some  plausible  theoretical 
noise  models  in  f-k  space  were  computed.  These  models  are  generated  from 

z 

combinations  of  the  |  crosspowers  |  which  results  from  disk  or  annulus 
distribution  of  power  in  the  K  plane  (at  a  fixed  frequency). 
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Channel  Coherences  for  25  March  1966  Noise  Sample 
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The  formulas  for  the  disk  and  annulus  are  as  follows: 


Disk: 


$  (f\  =  - j  (Zrrf  |x|) 

*12  TTf  X  J1  V 


where 


V  =  lower  velocity  or  radius  of  the  disk  (km/sec) 


f  =  frequency  (cps) 

X  =  separation  of  seismometers  (km) 


Annulus: 


*12<f> 


TTf  (X 


V2 

U‘ 


2TTf  lx 


_L  t  (ilLIaI 
v  M  v 

u  \  U 


where 


=  velocity  of  outside  radius  of  annulus 


=  velocity  of  inside  radius  of  annulus 


,  st 


Jj  =  1  -order  Bessel  function 


For  shifted  disk,  the  crosspow'er  from  a  centered  disk  is 
modified  by  multiplying  the  formula  for  the  centered  disk  by 

1 


where 


exp  i  2tt  ^  1 5? |  cos  ij/j 


is  the  vector  location  (in  the  K  plane)  of 
the  center  of  the  disk 


X 

* 


is  the  vector  separation  of  the  seismometers 

— *  — * 

is  the  angle  between  K  and  X 
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Figure  IV-4(a,b,  c)  shows  the  plots  calculated  for  three 
different  disks  in  the  K  plane.  The  model  used  was  a  disk  of  uniform 
power  density  out  to  wavenumbers  correSi  onding  to  the  velocities 
shown.  The  horizontal  scale  is  f  |X|,  allowing  the  frequency  to  be  scaled 
according  to  the  sensor  separation  desired.  Thus,  for  a  10-km  separation, 
the  horizontal  scale  would  be  interpreted  as  ranging  from  0  to  2.  6  cps  in 
frequency. 

Figure  IV-4  shows  the  coherence  obtained  from  an  annular 
noise  model.  The  noise  was  assumed  to  be  uniformly  distributed  in  an 
annulus  between  V  ~  10  km/sec  and  V  =  20  km/sec. 

Figures  lV-4e  and  IV-4f  show  the  coherences  obtained  from 
a  model  made  up  of  two  and  three  shifted  disks,  respectively.  These  disks 
were  designed  to  correspond  to  a  circular  region  in  the  K  plane,  which 
would  fit  precisely  into  an  annulus  of  12  km/sec  to  18  km/sec.  This  is 
equivalent  (in  size)  to  a  disk  centered  at  K  =  0  (V  =  ®)  with  a  minimum 
velocity  of  72  km/sec.  The  center  of  these  shifted  disks  is  a  V  =  14.4  km/ sec. 

For  models  involving  shifted  disks,  coherences  are  a  function 
of  the  relative  orientation  of  the  noise  power  model  and  the  sensor  locations. 
For  the  2-disk  model,  the  disks  and  the  sensor  are  centered  on  the  X  axis. 

For  the  3-disk  model,  the  disks  are  separated  by  120°,  with  one  disk  centered 
on  the  K  axis  and  the  sensors  placed  on  the  X  axis. 

A 

The  models  are  all  based  on  various  concepts  of  the  distribution 
of  the  mantle  P-wave  noise. 

The  two  larger  disks  (V  =  8  km/sec  and  V  =  12  km/sec) 

0  6 

treat  the  P-wave  energy  as  isotropic.  For  large  sensor  separations 
( I X |  >  10  km/sec),  these  models  would  give  essentially  no  coherence  except 
at  low  frequencies  (f  <  0.5  cps). 
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The  10-  to  20-km/3ec  annulus  model  approximately  represents 
a  model  of  P-wave  noise  which  is  uniformly  generated  from  teleseismic  dis¬ 
tances  to  the  core  shadow  region.  Again,  there  is  little  coherence  above 
0.  3  cps  for  |x|  >  10  km. 

The  disk  =  72  km/sec  represents  the  P-wave  noise  as 
almost  a  point  source.  This  would  be  a  good  model  if  the  P-wave  energy 
were  being  generated  by  a  single  storm  at  teleseismic  distances.  This 
model  gives  considerable  coherence  over  a  broad  band,  even  for  10-  to 
20-km  spacings. 

Figure  IV-4c  shows  the  results  of  two  and  three  point-like 
sources.  For  these  models,  the  interference  patterns  of  the  sources  cause 
the  coherences  to  have  peaks  and  nulls;  however,  the  peak  coherences  are 
still  significant  over  a  good  range  of  frequencies,  even  for  10-  to  20-km 
spacings. 

The  2-channel  coherences  (Figure  IV- 1)  are  essentially  0, 
even  for  spacings  near  10  km  (B  ring),  except  at  the  microseism  peak. 

The  coherence  near  the  microseism  pieak  is  likely  caused  by  the  surface- 
mode  noise  rather  than  the  P-wave  noise.  Therefore,  the  measured  2- 
channel  coherences  seem  to  preclude  that  the  P-wave  noise  is  due  essentially 
to  a  few  point-like  sources  unless  the  subarray  processing  is  destroying  the 
coherence.  Even  with  noise  dominated  by  point-like  sources,  little  coherence 
would  be  expected  at  large  distances  (50  to  100  km). 
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SECTION  V 

MULTICHANNEL  COHERENCE  OF  LASA  NOISE 

Multichannel  coherence  of  subarray  outputs  was  studied  for 
the  noise  samples  given  in  Table  r-l.  To  calculate  multiple  coherences 
between  AO  and  subarrays  through  the  D  ring  for  these  noise  samples,  the 
time  traces  were  filled  out  to  4096  points  by  adding  0's  and  estimating  the 
crosspower  matrix  from  the  direct  transform  method  using  the  Cooley- 
Tukey  algorithm.  In  all  cases,  the  smoothing  was  0.2  cps,  or  about  8C 
frequency  points  (Af  =  1  / (nAt)  =  1/409.6  cps). 


Table  V-l 
NOISE  SAMPLES 


Date 

Approximate 
Starting  Time 

No.  Points  Used 
(At  =0.1  sec) 

29  October  1965 

21:01 

3400 

4  November  1965 

00:42 

2950 

10  November  1965 

02:00 

3300 

25  November  1965 

03:40 

2048 

4  December  1965 

03:07 

3400 

21  December  1965 

08:41 

3490 

21  January  1966 

06:57 

4096 

5  February  1966 

03:02 

4096 

L5  March  1966 

04:34 

4096 

8  April  1966 

05:02 

4096 

1 5  pril  1966 

06:44 

4096 

29  April  1966 

09:34 

4096 
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As  indicated  in  Section  II,  such  smoothing  can  severely  reduce 
the  coherence  when  the  separation  of  sensors  is  fairly  large  (D  ring  or  greater). 
This  effect  should  be  reduced  somewhat  by  considering  the  multiple  coherences 
by  rings.  Some  sensors  on  the  D  ring  should  be  located  in  such  a  way  that 
they  have  relatively  little  movcout  with  respect  to  AO. 

Also  calculated  were  the  multiple  correlations  from  several 
reasonable  models  in  the  f-k  space  of  the  mantle  P-wave  noise.  These  are 
presented  for  their  intrinsic  interest  and  as  a  basis  for  interpreting  the  •' 

observed  multichannel  coherences. 

The  noise  sample  recorded  on  25  March  1966  served  as  the 
pilot  test  for  the  multichannel  coherence  studies  and,  therefore,  was  pro¬ 
cessed  more  intensively.  Using  this  4800-point  8-min  noise  sample,  the 
following  multichannel  filter  systems  were  designed: 

•  A  filter  to  predict  AO  MCF  output  from 
Bl,  B2,  B3,  C2,  C3,  C4,  Dl,  D2,  D3, 
and  D4  MCF  outputs 

•  A  filter  to  predict  AO  E  output  from  Bl, 

B2,  B3 ,  C2,  C3,  C4,  Dl,  D2,  D3,  and 
D4E  outputs 

•  A  filter  to  predict  AO  Z-21  output  from 
Bl,  B2,  B3,  C2,  C3,  C4,  Dl,  D2,  D3, 
and  D4  Z-21  outputs 

•  A  filter  to  predict  AO  MCF  output  from 
C2,  C3,  C4,  Dl,  D2,  D3,  and  D4  MCF 
outputs 

•  A  filter  to  predict  AO  MCF  output  from 
Dl,  D2,  D3,  and  D4  MCF  outputs 

All  of  the  above  filters  were  designed  in  the  time  domain 
from  prewhitened  data  and  were  41  points  (4.  1  sec)  long. 
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Figure  V-l  shows  the  fractional  prediction  error  as  a  function 
of  frequency  in  predicting  the  AO  output  using  subarrays  Bl,  B2,  B3,  C2,  C3, 
C4,  Dl,  D2,  D3,  and  D4  where  -3  db  corresponds  to  50-percent  predictability 
and  -6  db  corresponds  to  75-percent  predictability.  These  curves  are  calcu¬ 
lated  by  taking  the  ratio  of  the  power  spectrum  of  the  prediction  error  to  the 
power  spectrum  of  the  reference  trace. 
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Figure  V-l.  Prediction  Error  for  25  March  1966  Noise 
Sample  (B,  C,  D  Rings) 
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This  noise  is  most  predictable  near  the  microseism  peak 
and  appears  to  be  moderately  predictable  below  0.5  cps.  Between  0.5  and 
1.  25  cps,  the  noise  is  slightly  predictable;  from  0  to  about  30  percent  of 
the  power  is  predictable.  Above  1.25  cps.  the  noise  is  generally  unpredictable. 

Summation  and  MCF  outputs  give  comparable  predictability, 
while  single  seismometer  outputs  are  somewhat  less  predictable. 
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Superimposed  on  Figure  V-l  (bottom  graph)  is  the  prediction 
error  given  by  the  multiple  coherence;  4096  points  are  used  and  the  smoothing 
interval  is  0.25  cps,  giving  approximately  100  points  in  the  frequency  domain. 
There  is  general  agreement  between  the  two  different  methods. 


Figure  V-2  shows  the  fractional  prediction  error  obtained 
using  the  MCF  outputs  of  subarrays  composed  of  the  C  and  D  rings  and  of 
the  D  ring  only  to  predict  the  AO  output.  Results  indicate  that  the  predictability 
comes  from  the  B  and  C  rings. 
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Figure  V-2.  Prediction  Error  for  25  March  1966 
Noise  Sample  (C  and  D  Rings) 
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Figures  V-3  through  V-14  show  the  calculated  multiple  coher¬ 
ence  of  subarray  MCF  outputs  for  the  12  noise  samples  listed  in  Table  V-l. 
These  figures  show  the  prediction  error  (E)  as  a  function  of  frequency  when 
using  different  combinations  to  predict  the  AO  output.  The  subarray  outputs 
used  to  predict  AO  in  each  plot  are  identified  above  each  plot.  In  general* 
the  noise  shows  little  coherence  outside  the  B  ring.  Coherences  between 
AO  and  the  C  and  D  rings  are  generally  almost  0  or  small;  however,  the 
15  April  1966  noise  sample  (Figure  V-13)  is  an  exception. 

Maximum  coherences  between  AO  and  tne  subarrays  included 
in  the  B  ring  are  usually  about  0.  75;  75  percent  of  the  power  in  AO  output 
is  predictable.  The  maximum  coherence  occurs  near  the  microseism  peak 
and  is  very  likely  influenced  by  surface-mode  noise  which  is  'getting 
through"  the  subarray  velocity  filtering. 

In  general,  it  appears  that  coherence  attributable  to  the 
mantle  P-wave  noise  exists  only  within  the  3  ring  and  at  frequencies  below 
0=8  cps-  This  coherence  is  generally  0.25  to  0.70  in  this  interval. 

Several  noise  samples  show  coherent  peaks  at  higher  fre¬ 
quencies  (1.  5  <  f  <  2.0  cps),  which  are  very  likely  the  result  of  locally 
generated  noise. 

The  15  April  1966  noise  sample  (Figure  V-13)  appears  to  be 

anomalous.  While  the  low  frequency  (f  <  0.8  cps)  portion  of  this  sample 

appears  to  have  characteristics  similar  to  the  other  noise  samples,  a  quite 

coherent  peak  occurs  at  approximately  1.0  to  1.2  cps.  The  15  April  1966 

noise  sample  is  unique  in  that  a  point-like  source  of  high-velocity  noise  is 

present  near  1.  0  cps.  This  source  constitutes  about  80  percent  of  the  total 

power.  Characteristics  of  this  noise  sample  are  discussed  in  more  detail 
4 

in  another  report. 
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Multiple  Coherences,  21  December  1965 
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Figure  V-14.  Multiple  Coherences 


Figures  V-15  through  V-26  show  the  multiple  coherences 
calculated  from  some  models  of  the  high-velocity  noise.  The  theoretical 
crosspower  matricies  are  calculated  in  the  same  manner  as  those  in 
Section  IV.  Three  separate  models  are  considered. 

Figures  V-15  through  V-18  show  the  multiple  coherences 
from  a  uniform-disk  model  (in  the  K  plane)  with  a  minimum  velocity  of 
9  km/sec  and  uniform  power  distribution.  This  model  is  an  approximation 
of  uniformly  distributed  P-wave  noise.  The  multiple  coherences  are  shown 

withO-,  10-,  20-,  and  50-percent  random  noise  added  to  the  crosspower  matrix. 

The  plots  include  additional  subarrays  in  the  prediction  of  AO  in  the  sequence 
D4,  D3,  D2.  Dl.  C4.  C3,  C2,  Cl.  B4.  B3.  B2.  Bl;i.e.,  the  first  plot  uses 
D4  to  predict  AO,  the  second  uses  D4  and  D3  to  predict  AO,  etc. 

Figures  V-19  through  V-22  show  the  multiple  coherences 
calculated  from  a  model  (in  the  K  plane)  which  includes  a  uniform  disk  with 
a  9-km/sec  edge  velocity  upon  which  three  disks  with  radii  corresponding  to 
a  72-km/sec  edge  velocity  are  superimposed.  The  small  disks  are  centered 
at  36  km/sec  on  120°  radial  lines.  Ee.ch  small  disk  contains  the  same  total 
power  as  the  large  disk.  The  multiple  coherences  are  shown  with  0-,  10-, 
20-,  and  50-percent  random  noise  added  to  the  crosspower  matrix.  The 
plots  include  additional  subarray  outputs  in  the  prediction  of  AC  in  the  same 
sequence  as  described. 

This  model  is  an  approximation  of  a  P-wave  noise  field  which 
has  a  uniform  component  and  three  strong  point-like  source  components. 

Figures  V-23  through  V-26  show  the  multiple  coherences 
calculated  from  a  model  (in  the  K  plane)  which  includes  a  uniform  disk  with 
9-km/sec  edge  velocity  upon  which  th-ee  disks  with  radii  corresponding  to 
a  72-km/sec  edge  velocity  are  superimposed.  The  small  disks  are  centered 
at  36  km/sec  on  120  radial  lines.  Each  small  disk  contains  four  times  the 
total  power  of  the  large  disk. 
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Multiple  coherences  are  shown  with  0-,  10-,  20-,  and 
50-percent  random  noise  added  to  the  crosspower  matrix.  The  plots  include 
additional  subarray  outfits  in  the  prediction  of  AO  in  the  same  sequence  as 
described. 


This  model  approximates  a  P-wave  noise  field  which  has  a 
uniform  component  and  three  very  strong  point-like  source  components. 

Study  of  these  theoretical  multiple  coherences  and  comparison 
with  those  calculated  for  the  I  ASA  noise  suggest  the  following: 

•  Anisotropic  P-wave  noise  field  (all  lower-velocity 
energy  having  been  removed  by  the  subarray)  would 
give  strong  coherence  only  at  low  frequencies  (f  < 

0. 4  cps)  and  only  within  the  B  ring  of  subarrays 

•  Three  strong  point-like  sources  superimposed 
on  an  isotropic  P-wave  noise  field  would  be 
expected  to  give  useful  coherences  only  within 

/  the  C  ring  of  subarrays;  it  is  assumed  again 

that  the  subarrays  eliminate  the  low- velocity 
noise  and  do  not  distort  the  high-velocity  noise 

•  Three  very  strong  point-like  sources  superimposed 
on  an  isotropic  P-wave  noise  field  will  give  good 
coherence  between  the  AO  subarray  and  D-ring 
subarrays  when  the  low-velocity  noise  has  been 
effectively  suppressed  by  the  subarrays 

•  No  broadband  strong  point-like  sources  show  up  in 
the  P-wave  noise  as  seen  by  using  the  L.ASA 

•  A  strong  component  (at  ’east  20  percent)  of 
uncorrelated  noise  is  in  the  subarray  outputs 

This  incoherent  component,  it  should  be  noted,  may  be  partially 

5 

due  to  1.6-km/iec  isotropic  surface-mode  noise.  The  MCF  system  would 
pass  a  considerable  amount  of  this  energy.^ 

A  visual  fit  of  the  models  available  to  the  LASA  noise  coherences 
shows  that  the  one  giving  the  best  fit  appears  to  be  the  9-km/sec  disk  to  which 
has  been  added  20-  to  50-percent  random  noise. 
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SECTION  VI 


GROUP- COHERENCE  PROGRAM 


Given  two  sets  of  traces  (seismomete-  outputs),  it  is  desired 
to  construct  for  each  set  a  multichamel  filter  system  which  will  maximize 


the  2 -channel  coherence  between  the  two  MCF  outputs. 


The  problem  can  be  formulated  in  a  similar  fashion  for  the 
time  domain,  replacing  coherence  by  correlation.  Computationally,  the 
pr  jblem  is  much  more  tractable  in  the  frequency  domain. 


Development  of  the  group- coherence  theory  given  here  is  not 
mathematically  rigorous,  since  it  uses  the  concept  of  the  Fourier  transform 
of  a  stochastic  time  series  which  is  not  defined.  A  careful  mathematical 


development  could  be  given,  following  very  closely  the  development  of  the 
multiple  coherence  given  by  Koopmans. 


By  taking  the  direct  transform  of  N  channels  of  time-series 
data,  one  generates  (at  a  fixed  frequency)  N  complex  random  variables 
having  a  covariance  matrix  Q  ,  which  is  the  crosspower  matrix  If  the 
N- component  complex  random  vector  is  partitioned  into  two  sets 


(i) 


the  cross  power  matrix  is  partioned  into 
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It  is  desired  to  compute  linear  filters  which  will  maximize  the 
coherence  between  the  outputs  of  the  two  partitioned  sets.  If  a  is  the  vector 
of  complex  filter  weights  (at  the  fixed  frequency)  for  the  first  set  (X  )  and  F 
is  the  filter-weight  vector  for  the  second  set,  the  quantity  to  be  maximized 
can  be  written  as 


E  (a*  Xl)  (r  *  X2) 


E  [ (a  *  X 1 )  (a.  *  X " )  ]  E  [<r  *  X2)  (r  *  X2)  ] 


a*n17r  |' 


(a*  a)(r*n22D 


E  denotes  expected  value,  and  *  denotes  conjugate  transpose. 

The  problem  also  can  be  defined  by  maximizing  la^n^T  | 

subject  to  constraints  a  Q^cc  =  1  and  T  =  ®*nce  a  ^12  ”  *s 

complex  in  general,  the  inethod  is  to  maximize  the  real  and  imaginary  parts 
and  show  that  they  lead  to  the  same  solution. 

Using  the  Lagrange  multiplier  method,  we  define 


y  =  i/2(a!!e  n12r  +  r*  n  :a)  -  (i/2)  x(a*nna- 1)  -  (i/2)  (.l(r*n22r  - 1) 


where  (a*  ni2  T  +  T  *  fi^a)  is  twice  the  real  part  of  a*  * 


Lemma  1 


The  following  two  lemmas  are  necessary  to  solve  the  problem 


aaj  a  nn  a  =  2*Aal  ‘  Ba2* 


—  a  *  Cl,,  a  =  2  (Aa,  +  Ba.) 
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where 


nn  *  A  +  Bi  and  a  *  a1  +  a2i 


If  we  let 


-J-  +  iJL 

aal  *a2 


3a  (a*nila)  =  2  nna 


Lemma  2 


note  that 


■h  (“*n12r  +  r*n21a)  =  2n12r 


In  defining 


a  ,  .  _a_ 

9“l  ‘  Sa2 


0  =  -Jjl  ,  JLL  „  4x  ,  o 

aa1  aa2  3a 


Taking  the  derivatives  of  y  with  respect  to  a  and  T  and  setting 
them  equal  to  0,  we  arrive  at 


^12  ^  a  =  0 


ni2  a  -  n  n22  r  =  0 


(6-1) 


(6-2) 
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Similarly,  maximizing  with  respect  to  the  imaginary  part  of 
the  coherent;,  we  define 


/  =  fa*  n12  r  -  r  *  n21a)  -  x'  (a*  nna- 1)  -y  n'(r*n22r- 1) 


Using  the  lemma 


(a«n12r-  r* 


°21  °^j  '  1  n] 


and  setting  the  derivatives  of  y'  with  respect  to  a  and  T  equal  to  0,  we  derive 


in12  r  -  x'  n  a  =  o 


(6-3) 


Ql2  a  “  ^  n22  ^  =  0 


(6-4) 


Multiplying  Equation  (6-1)  by  a*  and  Equation  (6-2)  by  T*  and 
applying  the  constraints  a  a=  1  and  F*  fi22  T  =  1,  we  derive 


a*  01o  r  =  X 


r  n21  a  =  n 


In  other  words,  X  =  (i  and  the  coherence 


|a*ni2r|2 

<a*nua)(r*a,2r) 


=  n 
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li 

0 


D 

D 


D 

D 


or 


Multiplying  Equation  (6-2)  by  Cl  ^  *  we  get 


n22  Q2ia  ~  X  T 


n22  Q21a 


(6-5) 


Multiplying  Equation  (6-1)  by  X  ,  we  have 


x*  r  =  x  x '  n, .  a 


12 


11 


(6-6) 


Using  Equation  (6-5)  in  Equation  (6-6),  we  derive 


ni2  ^22  ^2!  a  =  1  X  fill  a 


which  is  a  generalized  eigenvalue  problem  that  can  be  solved  for  the  co¬ 
herence  (  |X  |2)  and  the  filtering  weights  (a),  which  are  the  components 
of  the  associated  eigenvector. 

Operating  similarly  on  Equations  (6-3)  and  (6-4),  we  arrive 
at 

-ia*  r  =  x' 


and 


12 


ir *  n  a  =  u 


Multiplying  Equation  (6-4)  from  the  left  by  Q  '  ,  we  have 

4  4 


in22_1  n21a 


and  multiplying  Equation  (6  —  3)  by  X  gives 


-i  *  n12r  *  x'  *  nna 


Using  Equation  (6-7)  in  Equation  (6-8),  we  derive 


^12  ^22  ^21  a  =  1^-  I  n  ii  a 


(6-7) 


(6-8) 


(6-9) 


which  is  the  same  problem  arrived  at  when  setting  the  derivatives  of  y  =  0, 
Thus,  maximizing  with  respect  to  the  real  or  imaginary  part  of  the  coherence 
is  equivalent  to  maximizing  with  respect  to  the  absolute  value. 


The  eigenvalue  problem  to  be  solved  is  thus  reduced  to 


^12  ^22  ^21  ’  ^  I  ^i  i^a  =  0 


(6-10) 


(Q11_1  ni2  n22_1  °21  '  lX|2jt)  a  =  °  (6-  i !) 

Note  that,  since  the  partition  of  a  covariance  matrix  is 


nll 

n!2 

n21 

^22 
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n  1 1 .  *,  n22,  and  n22  1  are  positive  definite;  also  since  ni2  *  fl21  *, 

(ni2n22-1  n2l)  is  positive  definite.  Since  (fl  ^  1  fl2  ,)  and  fl  ^  are 

positive  definite,  there  exists  a  nonsingular  matrix  F  such  that 


F  ^12  ^22  ^21^  F  = 


lx/  ...  0 


I* 


P'J 


=  A 


where  A  is  the  matrix  of  eigenvalues  of  Equation  (6-10)  arranged  so  that 
| X.  | £  |  X.  |  for  j  <  i  and 

F*nuF  =  I 

where  I  is  the  identity  matrix. 

Since  F*  C 1^  fl22  1  F  =  A  and  F  *  fl^  F  =  I,  we  have 


F*ni2n22‘ln21F=  F*nUFA 


and  thus 


ni2  n22  n21  F  =  nil  F^ 


(6 


Since  A  is  a  diagonal  matrix, 


n.,  n.  ~l  n_.  f.  =  Q, ,  f.  |x.  |2 

12  22  21  i  11  x  l 


where  f^  is  the  i**1  column  of  F  and  |  X^  |2  is  the  i^  eigenvalue. 


-12) 
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Thus,  F  is  the  matrix  with  columns  that  are  the  eigenvectors 
of  Equation  (6-10).  Note  that,  from  Equation  (6-12),  we  have 


nil  ^12  Q22  °21  =  F  A  F  1 


12  "22  "21 


(6-13) 


which  will  be  useful  in  the  following  equations. 

The  problem  [Equation (6- 10)]is  solved  by  an  iteration  commonly 
known  as  the  power  method.  Taking  a  starting  solution  vector  of 


a<0>  -(i . if 


and  defining 


a  (N  +  .1 )  = 


nil‘lni2n22'ln21a(N> 


[a*(N)a(N)] 


nil"‘  ni2n22‘‘n21  6(N) 


where 


8  (N)  =  —-a. (N) 

|a(N)| 


(6-14) 


the  result  of  N  iterations  is 


BIN)  =  cN(nH1  ni2  n22‘‘  n2,)Na(0) 


where 


N-l 

IT  |a, 

i=0  1 


=  Cn(fAf_1)N  a  (0) 


[from  Equation  (6-13)] 


=  cNf(A)NF-1a(0)=  CjH  I^F— ! L_  (A)NF-Ia(0| 
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As  N  -*  <*>,  this  approaches  the  matrix 


Thus,  we  obtain 

lim  CN  |Xj  |2N  F  p  0 
N  -*  * 

0  . 


CN  |X.|2N  a«0) 

n-cp  n1  r  i  (gj) 


n"  CN  U1|2N  ,£1)K 


(6-15) 


where 


fj  is  the  first  column  of  F 

~T  .  -1 

.  .  is  the  first  row  of  F 

~T 


K  is  set  equal  to  ^  ^  a(0)  in  the  second  equality 


From  Equation  (6-14), 

_  jjc 

1  =  B*(N)3(N)  =  CN  [f  AN  F"  1  a  (0)]  Cn[fAN  F‘‘a(0)] 

By  factoring  out  |\.  |  ~nd  taking  the  limit  as  N-»«,  we  have,  from 
Equation  (6-15), 

1 '  [suj-f  [<K,<V]*  [<K)«v] 

Thus, 


1 


lim 
N  -•  ® 


W 


2N 


K  |  |f 


and  6  (N)  approaches  a  constant  times  the  first  column  of  F  (the  largest 
eigenvector). 

The  program  actually  checks  the  magnitude  of  the  difference 
between  each  successive  iteration  of  the  eigenvector  and  stops  the  iterations 
when  the  magnitude  becomes  sufficiently  small.  Then,  the  final  iteration  is 
normalized  so  th?i  a  *  nn  a  =  1,  and  the  associated  eigenvalue  is  calculated 
from 

Q1 1  ni2  ^22  ^Zl  a  =  ^1  ^  a 
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The  T  filter  vector,  associated  with  the  second  group  of 
sensors,  is  calculated  from  Equation  '6-2)  by 

r  =  >-  n22  °‘21a 


where  X  is  the  square  root  of  |x|2.  It  should  be  noted  that  an  arbitrary 
factor  of  e^  is  associated  with  the  computation  of  the  filter  vector  T.  Since 
we  maximize  with  respect  to  |  \  |  ,  this  arbitrary  phasor  is  to  be  expected, 
for  it  does  not  change  the  coherence.  This  may  be  considered  physically  as 
being  equivalent  to  moving  the  output  point  of  the  MCF  system  designed  to 
maximise  the  coherence. 

This  procedure  can  be  continued  to  find  the  successively  smaller 
eigenvalues  and  their  associated  eigenvectors. 


^  1 1  ^12  ^22  n21  "  F^F 


-1 


has  the  form 


I>A?T 

i=i 


where  f.  is  the  ith  eigenvector  (the  il“  column  of  F)  and  f,*  is  the  l  row 


.th 


.th 


of  F  . 


Since  F*  0  F  s  l,  it  follows  that  F  1  =  F*  0^  or 


f.T  =  f.*n.  Thus, 

i  l  11 


hf1(Wnz.-lhl2  ■E'.'h'T 

i=2 


0  .  .  .  .  0 
.  lx,  I2 


0  |x 


'  .2 


-1 
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i  i2 

The  largest  eigenvector  of  the  modified  matrix  is  |  |  •  The  iterative 

procedure  can  again  be  applied  to  find  the  second  largest  eigenvalue  and 
associated  eigenvector. 

In  actual  calculations,  the  crosspower  matrix  0  is  estimated 
by  direct-transtorm  methods  and  smoothed. 

There  has  been  no  attempt  to  derive  the  distribution  of  the 
estimates  of  the  maximum  coherences  (the  eigenvalues).  However,  we  have 
made  some  estimates  of  the  measured  group  coherence,  under  the  hypothesis 
that  the  traces  were  uncorrelated,  by  running  the  program  on  traces  coming 
from  a  random  noise  generator;  the  results  are  summarized  in  Table  VI- 1. 
The  table  shows  that  a  large  amount  of  data  is  necessary  if  one  wants  to 
maintain  a  reasonable  smoothing  interval  and  have  reasonable  reliability 
(stability)  in  the  estimate  of  the  group  coherence.  Apparently  a  good  rule- 
of  thumb  would  be  to  require  10  independent  frequency  data  vectors  for  each 
channel. 


Table  VI- 1 

GROUP  COHERENCES  FOR  RANDOM  NOISE 


Number  of  Trace* 
in  Each  Partition 
(1/2  number  of 
total  trace*) 

Time 

Interval 

(pointa) 

(aec) 

Smoothing  Interval 
in  Frequency 
Domain 
<cp.) 

Number  of  Dtaic 
Frequency  Intervala 

Number  of  Eatimatea 
of  Oroup 
Coherence 

Calculated  J  K  |^'a 

19 

4096  pt 

409.6  aec 

0.2 

81 

3 

0,64.0.59,  0.64 

10 

4096  pt 

409.6  aec 

0.2 

81 

3 

0.  36,  0.  335,  0.295 

10 

8192  pt 

819.2  ««c 

0.2 

162 

1 

0,  155 

10 

12268  pt 

1228.8  aec 

0.2 

243 

3 

0.  12,  0. 12,  0.135 
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SECTION  VII 


GROUP-COHERENCE  PROCESSING  OF  LASA  NOISE 


Does  the  noise  include  any  elements  which  are  strongly  co¬ 
herent  among  subarrays  ?  This  is  the  basic  question  to  be  answered  by  group- 
coherence  processing  on  the  LASA  noise. 

Used  for  processing  was  the  27  April  1966  noise  sample  which 
consisted  of  three  8-min  tapes  forming  a  24-min  noise  sample  with  a  gap  of 
only  a  few  tenths  of  a  second  between  tapes.  Based  on  results  of  the  study  of 
group-coherence  estimates  using  random  traces  (Table  VI-1),  the  amount  of 
data  waj  cons’dered  to  justify  use  of  about  10  channels  only  from  each  sub¬ 
array  while  maintaining  reasonable  resolution  (0.  2  cps). 

Chosen  from  subarrays  AO  and  B4  were  10  channels.  Sensor 
locations  are  shown  to  scale  in  Figure  VII- 1.  As  a  control,  10  channels  were 
selected  from  subarray  F4  and  the  group  coherences  were  calculated  for  AO- 
64  and  A0-F4. 

Crosspower  matrices  are  estimated  by  taking  the  direct  trans¬ 
form  of  4096  points  from  each  of  three  noise  tapes  (which  have  been  pre¬ 
whitened),  forming  the  crosspower  matrix  for  each,  and  smoothing  over  81 
frequency  increments  (0.2  cps).  Then,  the  three  crosspower  matrices  are 
stacked.  Additionally,  the  following  are  computed  for  comparison: 

•  A0-B4  group  coherence  using  only  one  cross¬ 
power  matrix  and  a  stack  of  two  of  the  cross¬ 
power  matrices  just  described 

•  2-chainel  coherence  from  the  MCF  output  of 
subarrays  A0  and  B4,  which  is  calculated  using 
4096  points  and  a  0.  2-cps  smoothing  interval 

•  Multiple  coherence  between  seismometer 
84,  subarray  B4,  and  seismometers  21, 

41,  61,  81,  22,  23,  25,  26,  27,  36,  56, 
and  76  irom  subarray  A0  (the  sensors  lo¬ 
cated  nearest  seismometer  84  of  B4),  which 
is  calculated  from  4096  points  of  unwhitened 
data  using  a  smoothing  interval  of  0. 1  ^s 

•cwvn cs  ••rvleta  division 


VII- 1 


Figure  VII-2  shows  results  for  the  preceding.  The  following 
are  the  salient  features  of  the  figure. 

•  Group  coherence  for  the  two  subarrays 
(which  are  only  5  km  apart,  center  to 
center)  is  not  extremely  high  except  at 
very  low  frequencies  (microseism  peaks) 

•  LASA  noise  is  uncorrelated  among  sub¬ 
arrays  above  approximately  1  cps 

•  IviCF  outputs  exhibit  considerably  less  co¬ 
herence 

•  No  coherence  is  indicated  between  AO  and 
F4  subarrays 

Group  coherence  is  very  high  only  at  0.  2  cps.  This  is  prob¬ 
ably  caused  by  a  narrowband  directional  surface  mode.  5  Coherence  in  the 
frequency  range  where  the  mantle  P-wave  noise  is  generally  dominant  (0.4 

to  1. 0  cps )  falls  in  the  range  of  0.  3  to  0.  7  and  is  roughly  as  coherent  as  the  noise 
from  a  single  point- like  source  (Section  IV). 

It  is  difficult  to  compare  exactly  the  group  coherence  with  the 
2-channel  coherences  cited  in  Section  IV,  for  there  is  some  question  as  to  what 
the  sensor  spacing  should  be  when  one  considers  group  coherence.  The  group- 
coherence  filters  could  shut  off  all  except  the  nearest  two  sensors  and  thus  be 
equivalent  to  a  2 -channel  coherence  but  much  more  closely  spaced  than  the  sub¬ 
array  centers.  An  examination  of  the  filter  responses  indicates  that  this 
did  not  happen.  However,  it  is  difficult  to  estimate  how  much  effect  sub- 
array-element  spacing  should  have  when  comparing  the  2-channel  coherences 
with  the  group  coherences. 

Wavenumber  responses  of  the  two  filter  systems  (for  AO  and 
B4)  have  been  computed,  but  no  useful  interpretation  of  them  can  be  made  so 
they  are  not  presented. 
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I-MEASURED  NOISE  1  SEGMENT  GROUP 

II-MEASURED  NOI SE  2  SEGMENTS  STACKED  COHERENCE  A„-B 
EQ-MEASURED  NOISE  3  SEGMENTS  STACKED  V-B, 


LASA  noise  is  apparently  incoherent  among  subarrays  above 
approximately  1  cps.  Group  coherence  above  1.0  cps  agrees  well  with  the 
values  computed  from  the  random  noise;  this  is  about  the  same  frequency 
where  the  noise  appears  incoherent  within  a  aubarray.  ^ 

The  MCF  processor  gives  an  output  which  shows  considerably 
less  2-channel  coherence  than  that  of  the  group-coherence  filters,  but  this  is 
to  be  expected.  The  group-coherence  filter  systems  should  "latch  on"  to 
elements  of  the  noise  field  which  are  most  similar  between  the  two  subarrays 
and  should  try  to  reject  the  remainder. 

Additionally,  the  group  coherence  is  greater  than  the  multiple 
coherence  between  the  outermost  element  of  the  B4  array  and  the  12  closest 
elements  from  subarr^y  AO.  This  suggests  that  there  is  some  coherence 
gain  by  rejection  of  part  of  th.  less  coherent  noise. 

Group  coherence  for  the  A0-F4  subarrays  stays  very  near  the 
level  expected  for  uncorrelated  traces.  There  could  be  some  coherence 
among  these  subarrays,  however,  which  would  r.vA  be  seen  due  to  the  smooth¬ 
ing  involved  in  the  estimate  of  the  crosspower  matrix.  These  subarrays  are 
so  widely  separated  that  even  a  perfectly  coherent  plane  wave  would  show  no 
coherence  unless  the  apparent  velocity  between  the  two  arrays  were  very  high 
(Section  II). 

The  A0-F4  group  coherence  ha^  significance  only  as  a  control 
to  indicate  that  the  calculation  of  the  group  coherence  is  behaving  in  a  reason¬ 
able  manner. 

A  large  signal  has  been  processed  as  a  check  of  group  coherence. 
Group  coherence  is  calculated  between  subarrays  B2  and  C2  for  a  signal  from 
Andreanof  Islands  (51.4°N,  179.  7°W,  A  =  46.5°,  Az  =  304.1°,  recorded  on 
23  November  1965).  From  the  LASA  sife,  this  would  give  a  A  of  46.  5°  and 
an  azimuth  of  304°.  The  apparent  horizontal  velocity  of  the  P  wave  would  be 
14.  3  km/  sec.  The  event  was  large  (mag  5.  6),  complex,  and  well- recorded  at 
LASA. 
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Group  coherence  is  calculated  using  13  sensors  from  each  sub¬ 
array  (10,  3  ring,  5th  ring,  7th  ring,  and  8th  ring).  Figure  VII-3  shows 
some  representative  traces  of  the  event.  Group  coherence  is  calculated  at 
1.  0  cps  using  3900  points  with  a  0.  2-cps  smoothing  interval.  To  stabilize 
the  inversion  problem,  5-percent  random  noise  is  added  to  the  crosspower 
matrix. 


Figure  VII-3.  Andreanof  Islands  Earthquake 


The  calculated  group  coherence  a;  1  cps  is  0.  96.  A  pair  of 
exactly  similar  traces  time- shifted  appropriately  and  processed  as  above 
(same  smoothing  and  random-noise  addition)  would  have  a  2-channel  coherence 
of  approximately  0.  82.  Recall  that  26  total  channels  have  been  used  to  pro¬ 
duce  the  output  which  gives  the  0.96  coherence  estimate;  this  means  that  the 
traces  can  be  used  to  reduce  random  noise.  Also,  there  are  approximately 
two  degrees  of  freedom  per  filter  point,  which  suggests  a  biased  high  estimate 
of  the  group  coherence. 


Figuve  VII-4  shows 
the  wavenumber  response  of 
the  group- coherence  filters 
for  subarrays  AO  and  B4.  The 
earthquake  signal's  location 
is  shown  by  a  cross  in  each 
plot.  These  plots  indicate 
that  the  group- coherence 
filter  systems  pr“?s  the  P- 
wave  energy,  which  would  be 
expected  to  be  the  most  co¬ 
herence  element.  The  group- 


coherence  filters  have  indi¬ 


cated,  in  this  simple  case, 
the  location  (in  f-k  space)  of 


the  most  coherent  element 


of  the  data. 


Figure  VII-4.  Responses  of  Group-Coherence 
MCF  Systems  at  1.  0  cps 
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A 


SUBAt’RAY  B2 


SECTION  VIII 
LASA  COHERENCE 


The  subarray  outputs  at  LASA  are,  in  general,  essentially 
incoherent  except  within  the  B  ring  of  subarrays;  even  here,  the  coherence 
is  generally  0.5  or  less  except  at  the  microseism  peak.  This  indicates  that 
there  can  be  no  important  gain  from  coherence  processing  of  subarray  outputs. 

Processing  designed  to  maximize  the  coherence  between  two 
aubarraye  9  km  apart  tends  to  indicate  that  the  P-wave  noise  contains  some 
moderately  coherent  elements.  This  experiment  used  insufficient  data  to 
be  truly  definitive,  however. 

The  obser'  jd  coherence  among  subarrays  could  be  strongly 

influenced  by 

•  The  geometiy  of  the  large  array.  (The 
filling  of  the  'unit  cell"  in  f-k  space  by  the 
seismic  energy  destroys  coherence) 

r,  Distortion  of  plane-wave  seismic  energy  by 
local  inhomogeneities  in  the  crust  or  upper 
mantle  (which  results  in  a  smearing- out  of 
energy  in  wavenumber  and  reduced  coherence) 

•  Subarray  processing  which  could  possibly  re¬ 
duce  P-wave  coherence 

•  Coherence  estimates  which  tend  to  be  biased 
on  the  low  side 

LASA  has  no  true  unit  cell  in  f-k  space  due  to  nonuniform 
spacing,  but  the  effect  is  qualititively  the  same  as  that  of  an  equally  spaced 
array.  The  model  studies  (Section  VI)  indicate  that  good,  fairly  broadband 
coherence  can  be  expected  out  through  the  D  ring  when  the  P-wave  noise  is 
strongly  dominated  by  point-like  sources.  If  the  P-wave  noise  were  isotropic, 
one  would  expect  to  find  essentially  no  coherence  above  the  microseism  peak 
for  the  LASA  spacing. 
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The  "random"  noise  component  has  been  observed  to  increase 

5 

(in  percentage  of  total)  with  frequency  within  subarray  dimensions.  This  is 
probably  the  result  of  the  smearing-out  of  plane  waves  by  local  inhomogeneities 
in  the  crust  or  upper  mantle.  These  effects  become  increasingly  severe  as 
the  wavelengths  decrease  (increasing  frequency). 

Such  differences  are  represented  on  a  larger  scale  by  the 
travel-time  anomalies  o  served  over  the  LASA  aperture.  These  large  travel¬ 
time  anomalies  strongly  suggest  that  the  coherence  of  the  noise  (unless  it 
were  a  point  source)  would  be  severely  reduced  by  local  inhomogeneities  in 
the  crust  or  upper  mantle  under  LASA. 

7 

There  is  strong  evidence  at  TFO  that  the  P-wave  noise  is 
(sometimes,  at  least)  strongly  dominated  by  somewhat  point-like  sources  such 
as  those  used  in  the  model  studies  (Sections  IV  and  V).  However,  no  coherences 
of  the  order  expected  from  such  a  model  were  observed  in  the  LASA  noise, 
again  suggesting  that  local  variations  in  crust  or  upper  mantle  affect  measured 
coherence. 

There  is  relatively  little  information  upon  which  to  assess  the 
effects  of  the  subarray  processing  on  coherence.  However,  for  the  one  case 
where  prediction  among  subarrays  was  made  with  a  single  seismometer,  with 
summation  of  subarray  elements,  and  with  an  isotropic  multichannel  filter 
system  (Section  V),  the  summation  and  MCF  outputs  gave  very  similar  results 
and  the  single  seismometers  gave  noticeably  poorer  prediction  results.  By 
no  means  is  this  a  definitive  experiment,  but  it  tends  to  suggest  that  the  sub¬ 
array  processing  does  not  severely  reduce  the  P-wave  coherence. 

The  group  coherence,  however,  was  considerably  higher  than 
the  comparable  2 -channel  coherence  for  the  27  April  noise  sample  (Section  V). 
One  possible  explanation  (of  several)  for  this  is  that  the  subarray  processing 
is  destroying  useful  P-wave  noise  coherence. 
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Data  available  to  date  are  insufficient  to  state  conclusively 
whether  subarray  processing  is  destroying  useful  P-wave  coherence. 

The  biasing  effect  of  smoothing  on  coherence  has  already  been 
discussed.  Coherence  estimation  is  made  difficult  by  the  degree  of  smoothing 
needed.  A  small  amount  of  smoothing  reduces  bias  and  prevents  the  over¬ 
lapping  of  noise  modes;  a  large  amount  of  smoothing  gives  reliability  to  the 
estimates.  This  has  been  a  particular  handicap  to  group-coherence  processing 
where  a  large  number  of  sensors  are  desirable  and  available  (25  from  each 
of  2  subarrays).  A  noise  sample  sufficiently  long  for  good  estimates  of  the 
most  coherent  elements  between  two  moderately  separated  subarrays  (e.  g. , 

AO  and  the  D  ring)  is  desirable.  This  would  require  approximately  10  hr  of 
data  if  the  analysis  were  to  include  surface -mode  energy,  or  about  4  hr  of 
data  if  one  were  interested  only  in  P-wave  noise. 

One  reason  for  lack  of  coherence  over  larger  distances  (lateral 
inhomogeneities  in  the  crust  and  upper  mantle)  may  be  a  function  of  location. 
The  Montana  LASA  is  located  on  a  thick  section  of  Cretaceous  and  Tertiary 

g 

shales.  Douze  has  suggested  that  ambient  noise  is  less  interpretable  at 
sites  where  thick  shale  sections  exist.  A  large  array  in  a  shield  area  could 
possibly  exhibit  better  P-wave  coherence  than  does  the  LASA. 
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13  ABSTRACT 

Much  of  the  analysis  of  the  mantle  P-wave  noise  proposed  for  this 
contract  depended  on  the  ability  to  perform  coherent  processing  by  treating 
seismic  noise  as  a  compilation  of  plane  waves  using  the  large  aperture 
available  at  LASA.  The  subarray  outputs  should  be  P-wave-limited,  and  the 
large  array  could  be  used  to  dissect  the  P-wave  noise;  therefore,  it  is 
extremely  important  to  understand  the  coherence  of  the  subarray  outputs. 

This  report  presents  the  results  of  coherence  studies  among  subarrays. 
Principal  results  of  the  study  include  the  development  of  a  program  to  compute 
multiple  coherences  and  thus  estimate  multichannel  coherence.  A  new 
processing  technique  to  design  filters  that  maximize  the  coherence  between 
the  output  of  two  groups  of  sensors  (group  coherence'  was  developed  and  pro¬ 
gramed.  Presented  are  results  of  measured  2-channel  and  multiple  coherence, 
indicating  that  the  subarray  outputs  generally  are  moderately  coherent  below 
about  1  cps  only  within  the  B  ring  of  the  subarrays.  Also  presented  are 
2-channel  and  multiple  coherences  for  several  possible  noise  models.  A  study 
of  these  coherences  suggests  that  the  most  reasonable  models  would  give  little 
coherence  between  AO  and  the  E  and  F  rings.  Comparison  of  model- study 
results  with  measured  coherences  suggests  that  a  considerable  portion 

l2^=pe?,Cent  °r  more)  of  the  subarray  output  power  is  not  interpretable  as 
seismic^ener^.  r 
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